Topological Description of the Aging Dynamics in Simple Glasses 
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We numerically investigate the aging dynamics of a monatomic Lennard- Jones glass, focusing on 
the topology of the potential energy landscape which, to this aim, has been partitioned in basins of 
attraction of stationary points (saddles and minima). The analysis of the stationary points visited 
during the aging dynamics shows the existence of two distinct regimes: i) at short times, t < t c , the 
system visits basins of saddles whose energies and orders decrease with t; ii) at long times, t > t c , 
the system mainly lies in basins pertaining to minima of slowly decreasing energy. The dynamics for 
t > t c can be represented by a simple random walk on a network of minima with a jump probability 
proportional to the inverse of the waiting time. 



PACS Numbers : 61.20.Lc, 64.70.Pf, 02.70.Ns 

The very slow approach to equilibrium of glassy sys- 
tems is a very interesting and intriguing physical pro- 
cess that, in the last few years, has received considerable 
attention from analytical, numerical and experimental 
point of view When a glassy system is suddenly 
brought in the glassy state from an high temperature 
(or low density) equilibrium state, it starts to explore 
the phase space with an off-equilibrium dynamics which 
slows down as the system "ages" . During this aging 
regime one time physical quantities, that are time in- 
dependent on average at equilibrium, slowly approach 
their equilibrium values, while two time quantities, like 
the correlation functions, depend on both times due to 
the lack of time translation invariance. The study of 
the off-equilibrium dynamics of glassy systems, and in 
particular of spin-glasses, has allowed to obtain impor- 
tant informations on the phase-space properties of the 
systems itself (behavior of the order parameter in the 
low temperature equilibrium phase |^-^1) and to formu- 
late off-equilibrium relations generalizing the usual equi- 
librium ones (generalized fluctuation-dissipation relation 
with the introduction of an effective temperature 
[3|. The investigation of the off-equilibrium properties 
seems to indicate the correctness of the conjecture stat- 
ing the similarity of structural glasses with some spin 
glass model (spin-glasses with one step replica symmetry 
breaking) The way a glassy system goes towards 

equilibrium is then not only an interesting problem by 
itself but it is also of great relevance in view of a deeper 
understanding of the nature of glass transition. In this re- 
spect, many efforts have been devoted to the understand- 
ing of the role played by the potential energy landscape 
in the equilibrium dynamics of supercooled liquids [ jl2| . 
The trajectory of the representative point in the 3iV con- 
figuration space can be mapped into a sequence of locally 
stable points (the so called Inherent Structures, IS p3|), 
that are the local minima of the total potential energy V: 
to each instantaneous configuration during the dynami- 
cal evolution of the system one can associate an IS by a 



steepest descent path in the V surface. The properties of 
the IS has been found to be very useful to clarify many 
features in the dynamics and the thermodynamics of su- 
percooled liquids in jL4]-[l7) and off Jl8|,[l9| equilibrium. 
An extension of such an analysis including all station- 
ary points of V (minima and saddles) |po|-p2| has been 
recently proposed to form a better representation of the 
system during the equilibrium dynamics. This approach 
gives a novel point of view on the equilibrium dynam- 
ics of supercooled liquids based on a detailed topological 
analysis of the potential energy landscape. In particu- 
lar it emerges [^lj that the order of the sampled saddles 
during the equilibrium dynamics is a well defined increas- 
ing function of temperature and it appears to vanish at 
the so-called mode-coupling p5[ | transition temperature 
Tmct, which then marks the crossover from a dynamics 
between basins of saddles (T > Tmct) to a dynamics 
between basins of minima (T < Tmct)- Similar results 
have been found p3| in the framework of Instantaneous 
Normal Modes analysis pi) . 

In this Letter we numerically investigate the off- 
equilibrium dynamics of a simple model glass, by map- 
ping the true dynamics in the configuration space on to 
the stationary points of V. In our glass the particles in- 
teract via a Lennard-Jones potential modified in such 
a way as to avoid crystallization [ p6| . The system is 
brought in an off-equilibrium initial condition by a sud- 
den isothermal density jump (crunch) across the liquid- 
glass transition line. We find a rather sharp cross-over 
between two dynamical regimes: i) a short times regime 
where the system point travels between basins of saddles 
of decreasing energy and order and ii) a long time regime 
where the system point moves between basins of saddles 
of zero order (minima) and decreasing energy, crossing 
basins of saddles of low order. The existence of two well 
separated dynamical regimes shows up also in the Mean 
Square Displacement (MSD). Here a logarithmic behav- 
ior is found at long times suggesting the picture of a sys- 
tem point performing a random walk between basins of 
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minima with a jump probability per unit time inversely 
proportional to time elapsed since crunch. This picture 
is quantitatively supported by the very good agreement 
found between the MSD measured along the MD tra- 
jectories and the MSD reconstructed via a random walk 
model. 

The system under study is a monatomic Modified- 
Lennard- Jones (MLJ) |2(| that is the usual 6-12 Lennard- 
Jones pair interaction (truncated and shifted at R c = 2.6, 
standard Lennard- Jones units are used hereafter) plus 
a small many body term which inhibits crystallization 
with negligible corrections to the equation of state of 
a monatomic Lennard- Jones system (see fll]j2(| for de- 
tails) . The system is composed of N = 256 particles en- 
closed in a cubic box with periodic boundary conditions. 
The initial liquid equilibrium configurations (po = 0.95, 
To = 0.5) are prepared by standard isothermal MD. 
At time t — the density is suddenly increased to the 
value p\ — 1.24, where the glass transition temperature 
(Tg(pi) ~ 1.4 (ll)) is well above the simulation tempera- 
ture To The subsequent evolution of the system is gen- 
erated up to time t = 3 ■ 10 4 by means of standard mi- 
crocanonical MD with time step r = 0.01. For each MD 
trajectory we quenched a set of logarithmically equally 
spaced configurations to the corresponding IS and saddle 
following the steepest descent path respectively on the 
potential energy surface V in the former case and on the 
auxiliary potential W = ||W| 2 (23-|| in the latter. 




FIG. 1. Energies as a function of time during off equilib- 
rium dynamics. Instantaneous configuration potential energy 
per particle subtracted by the constant term |To (open cir- 
cles), inherent structure saddle energy e s (full triangle) and 
the inherent structure minimum energy e m (open diamonds) . 
Error bars are shown only when larger than symbol size. The 
arrow indicates time t c where the crossover from dynamics be- 
tween basins of saddles to dynamics between basins of minima 
takes place. 

In Fig.|]we report in a logarithmic scale the time evolu- 
tion of: a) the instantaneous configuration potential en- 
ergy per particle subtracted of the constant term |To 
(open circles), b) the saddle's energy e s (full triangles), c) 



the IS energy e m (open diamonds) . The remarkable coin- 
cidence between e s and — §To suggests to think of the 
instantaneous energy ej as the sum of two contributions: 
the saddle's energy (slowly approaching equilibrium) and 
the energy of nearly N fast harmonic equilibrated degrees 
of freedom. It is evident the presence of a rather sharp 
transition in dynamical behavior of ej and e s at time 
t c ~ 20. In particular the time evolution of e, for t > t c 
seems to be driven by the inherent dynamics between 
minima as suggested by the underlying time evolution of 

■ 

In Fig.| the mean saddle's order n s (number of nega- 
tive eigenvalues of the potential energy Hessian computed 
on the saddle configurations) is reported as a function 
of time. The dynamic transition observed at t c in the 
energy time evolution in Fig.[l] corresponds here in the 
mean order n s falling below the value n s = 1 or, in other 
words, for times longer than t c the system point starts 
to explore extensively the basins of saddles of order zero 
(minima). The inset in Fig.|| shows the order n s com- 
puted along a single MD trajectory. One can see that for 
times t > t c the systems visits basins of minima passing 
through basins of low order saddles. 
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FIG. 2. Mean inherent saddle order n s (number of negative 
eigenvalues of the Hessian matrix of potential energy evalu- 
ated at the inherent structure saddle) as a function of time 
t after crunch. The inset shows n s computed along a single 
trajectory. 

In a recent study of the MLJ landscape probed in equi- 
librium pl[ some of us found that the mean elevation 
energy of a saddle of order n s from the underlying min- 
imum is proportional to n s (this is at variance with the 
result reported in [ p2[ , where a linear relation is observed 
between n s and the absolute energy of the saddles). This 
is clearly a property of the energy landscape itself and 
one expects that the same properties are also valid in the 
off-equilibrium dynamics. 

In Fig.|| where we report e s together with e m + 
n s AE/N with AE a fit parameter. This result confirms 
the picture of an energy landscape where the saddles of 
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order n s + 1 lay above the saddles of order n s of a con- 
stant quantity AE. From the data we obtain the value 
AE = 10. The value of AE could also be obtained 
from an equilibrium analysis, as in |2l| where a value 
of AE — 3.6 is found for the case p = 1. We have not 
carried out such an equilibrium analysis for the density 
p = 1.24 used here. 



-5.00 
-5.05 
-5.10 
-5.15 
-5.20 
-5.25 



T] 1 1 I I I I 1 1 1 1 1 I I I I I l| 1 1 I I I I ■ 1 1 1 1 I I I I 1 1 1 



e +n AE/N 

m s 



- □ 



10" 



10' 



10 



10" 



t 



FIG. 3. Saddle's energy e s (solid symbols) evaluated dur- 
ing off equilibrium dynamics compared with the energy of 
the minima e m increased by the quantity n s AE/N (open 
symbols) where n s is the time dependent saddle's order and 
AE — 10 is a fitting constant representing the mean energy 
separation between a saddle of order n s + 1 and a saddle of 
order n s . 
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FIG. 4. Mean square displacement r 2 (t,t w ) as a function 
of time t for four selected t m values: 1,5,20 and 200 (solid 
symbols). The straight lines are the prediction of a random 
walk model between minima with jump probability per unit 
time decreasing as the inverse of time, p(t) = c/t. According 
to the model the slope of these lines are d^c — 0.019 where 
do is the mean square distance between adjacent minima (see 
text below for details) . The arrow indicates time t c where the 
crossover from dynamics between basins of saddles to dynam- 
ics between basins of minima takes place. The system size is 
N = 256. 



Another interesting quantity in describing the dynami- 
cal processes is the mean square displacement of the con- 
figuration at time t from that at time t w after crunch: 



r 2 (*,««,) 



Ti(t w )\ 2 ) , 



(1) 



where (...) is an average over initial configurations. In 
Fig.|| we report r 2 (t,t w ) as a function of time t for four 
selected t w values: t w = 1,5,20,200. Also in this case 
t c ~ 20 marks a transition in the dynamical behavior. In 
particular the long time regime shows a well defined loga- 
rithmic dependence on time. This suggests the picture of 
the diffusion process proceeding through elementary un- 
corrected events occurring with a probability per unit 
time p(t) ocl/t [^7| . As previously observed, in the long 
time regime the system point travels between basins of 
minima. It is therefore natural to identify the elemen- 
tary diffusion events as jumps between adjacent basins 
of minima. In order to give further support to this con- 
jecture we calculated the cumulative number of jumps 
Nj umps {t) defined as the number of jump events in which 
the constantly monitored energy of the IS makes a sharp 
transition. 
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FIG. 5. a) Cumulative number of jumps between adjacent 
minima visited by the system during the off-equilibrium dy- 
namics; the solid line is the best fit for t > 20 with slope 
c = 1.0. b) Mean square distance between adjacent minima 
d 2 ; the solid line is the constant value dg ~ 0.019 reached for 
times greater than t c ~ 20. The system size is N = 256. 

The quantity Nj urnps is plotted in Fig.^(a). For times 
t > t c N jumps displays a nice logarithmic behavior that 
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is what one would expect if the probability per unit time 
of a jump event is pit) = c/t. In other words the time 
needed to escape from an IS reached by the system after 
a time t w from crunch is proportional to time t w itself (a 
similar behavior is found in the off-equilibrium dynam- 
ics of a ferromagnetic Ising system where the persistence 
time of a bubble in which spins have a similar value is 
proportional to time t w |2|,|2^]). In Fig||(a) the solid line 
represents the fitting function clog(i) with c = 1.0. In 
Fig.^(b) we report the mean square distance d 2 between 
adjacent minima probed by the system during the aging 
dynamics: 



d\t) 



(t + ) 



(t-)\ 



(2) 



where t is the time of the jump between two different 
basins of minima and rj s (t~), r* s (t + ) are the coordi- 



nates of IS respectively before and after time t. The 
quantity reported in figure approaches the long time con- 
stant value g?q ~ 0.019. For uncorrelated jumps, the 
mean square displacement would be simply given by 



r 2 (t,t w ) = dlclog{t/t v 



(3) 



This is in a remarkably good agreement with the MD 
data as shown in Fig.^J, where we draw straight lines with 
slope c?qC = 0.019 (as obtained from the values of d^ and 
c derived from Fig.||) and intercept chosen as to best fit 
the data for t > 20. 

In conclusion we have numerically studied the off- 
equilibrium dynamics of a monatomic Modified-Lennard- 
Jones system by sudden isothermal density jumps from 
the liquid to the glassy phase. Analyzing the time evo- 
lution of the energy and the mean square displacement 
we found two well defined dynamical regimes, which can 
be related to the different ways the system explores the 
potential energy landscape: i) at short times the system 
point moves downward passing through basins of sad- 
dles with decreasing energy and order, and ii) at longer 
time the relevant dynamical processes are jumps between 
basins of minima with slower decreasing energy. In the 
whole time regime the instantaneous potential energy 
can be thought as due to a sum of two contributions: the 
energy of the underlying saddle's configurations e s plus a 
nearly harmonic term of the fast degrees of freedom 3/2T. 
Moreover the linear energy gap between saddles and un- 
derlying minima e s — e m as a function of saddles order n s 
supports the picture of a simple energy landscape with 
only a single energy barrier parameter AE which repre- 
sents the energy gap between a saddle of order n s and 
a saddle of order n s + 1. The relevant processes in the 
long time aging regime are well represented by a simple 
random walk model with elementary events correspond- 
ing to jumps between minima with a probability per unit 
time proportional to the inverse of time elapsed since the 
system has been brought out of equilibrium. 
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